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Remarks on the realization of the Atiyah-Singer index theorem in lattice 
gauge theory* 

C. R. Gattringer^t^ i. Hip'' and C. B. Lang^ 

'^Department of Physics and Astronomy, University of British Columbia, Vancouver B.C., Canada 
'^Institut fiir Theoretische Physik, Universitat Graz, A-8010 Graz, Austria 

We discuss the interplay between topologically non-trivial gauge field configurations and the spectrum of the 
Wilson-Dirac operator in lattice gauge theory. Our analysis is based on analytic arguments and numerical results 
from a lattice simulation of QED2. 
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1. MOTIVATION 

It has been known ^ for more than two 
decades that topologically non-trivial configura- 
tions play an essential role for a proper under- 
standing of QCD. However, it is not straightfor- 
ward to go beyond a semiclassical analysis when 
implementing topological ideas in the continuum 
path integral. There are at least two reasons: 
Firstly the continuum path integral (if it can be 
constructed at all) has support only on configu- 
rations without the necessary smoothness condi- 
tion for a definition of the topological charge - in 
general the gauge fields contributing to the con- 
tinuum path integral are simply too rough to be 
classified with respect to their topological charge. 
Secondly one would also like to use the Atiyah- 
Singer index theorem ||] in the fully quantized 
theory. It relates the topological charge of a dif- 
ferentiable background configuration to the zero 
modes of the Dirac operator. Again the lack of 
smoothness of the gauge fields contributing to the 
path integral does not allow the application of the 
index theorem beyond a semiclassical analysis. 

On the lattice the situation is different. There 
are several definitions of the topological charge 
on the lattice which allow to classify (almost) all 
gauge field configurations contributing in the con- 
tinuum limit. In contrary to the continuum a 
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decomposition of the path integral into topolog- 
ical sectors is possible. Concerning the Atiyah- 
Singer index theorem, the situation is more diffi- 
cult. There is no analytic result for an index the- 
orem on the lattice, but it is generally believed 
that an equivalent result should be manifest on 
the lattice in a probabilistic way. 

Although several aspects of the interplay bet- 
ween topologically non-trivial configurations and 
the spectrum of the Dirac operator on the lat- 
tice have been analyzed due to the immense 
computational cost no complete study in fully 
quantized QCD4 has been accomplished. The 
numerical data which we present in this contri- 
bution come from a study [Q in lattice QED2. 
This model has similar structural features as 
QCD4 (topologically non-trivial configurations, 
existence of the anomaly, U(l)-problem etc.) but 
is computationally much less demanding. This 
allows for a proper understanding of the role of 
topology in that simple model, and indicates what 
to expect for a lattice study of QCD4. 



2. SYMMETRIES OF THE WILSON- 
DIRAC OPERATOR 

Before numerically analyzing the spectrum of 
the Wilson-Dirac operator M it is helpful to re- 
assess the symmetry properties of the operator 
(fermion matrix) and thus of the spectrum. We 
denote the fermion matrix as 7\f = 1 — kQ, where 
K is the hopping parameter and Q the hopping 
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matrix 

±D 

Q{x,y) ^ ^ i^+ll^Wiyix ~ {>) Sa:-i>,y , (1) 
u=±l 

where we use the notation 7_y = —j^, and 
U-v{x -\- v) = Ui,{x)'^ . D is the number of di- 
mensions (2 or 4) and Uu{x) are elements of the 
gauge group (U(l) for D = 2 and SU(N) for D 
— 4). We work on a finite L^~^ x T lattice with 
mixed periodic boundary conditions. The exis- 
tence of the following similarity transformations 
is well known 

rsOTs = Qt^ SQS = -g , (2) 

where the latter only holds for even L and T. The 
matrices Ts = 75^a;,y and S = (_i)a;i+..+2;D ^ 
are both unitary and hermitian. This implies 
that the spectrum of Q is symmetric with respect 
to reflection on both the real and the imaginary 
axis. Thus the eigenvalues of the hopping matrix 
on a lattice with even L and T come in complex 
quadruples or in real pairs. 

The special role of the real eigenvalues A of the 
hopping matrix (and thus the full fermion matrix 
M = 1 — kQ) can be appreciated only by under- 
standing the chiral properties of their eigenvec- 
tors. Denote by v\, the eigenvectors of Q with 
eigenvalues \,^. Using the symmetries (^) it can 
be shown Q that 

v^^^T^vx ^0, (3) 

only for /i = A. In particular for the diagonal 
elements (this result can already be found in |^) 

WA^ Ts «A ^ , (4) 

only for real eigenvalues A. This result has to 
be compared to the continuum result {tp denotes 
some eigenstate of the continuum Dirac operator) 

(V',75V') ^0, (5) 

only if is a zero mode. This indicates that 
only the eigenvectors of the fermion matrix cor- 
responding to real eigenvalues can play the role 
of the zero modes of the continuum. 

At this point we would like to remark that for 
an analysis of the realizations of index theorems 



on the lattice the Wilson form of the lattice Dirac 
operator is superior to the staggered version. The 
latter is antihermitian and thus has only purely 
complex eigenvalues lacking a straightforward cri- 
terion for identifying zero modes. Furthermore it 
is also not clear how to relate the (real) eigenval- 
ues of the hermitian operator T^Q to the spec- 
trum of Q with its more sophisticated features. 

3. NUMERICAL RESULTS IN QED2 

Having established the real eigenvalues of the 
fermion matrix as the trace of gauge field configu- 
rations with non-trivial topology one can simplify 
a numerical analysis by concentrating on the real 
eigenvalues alone. Here we discuss results from 
a numerical simulation of QED2 on the lattice 
Q. In this model the evaluation of the geomet- 
ric definition |^ of the topological charge v\U] is 
straightforward, giving 

^[^] = ^E^^W e (6) 

The plaquette angle Op{x) is introduced as 
Op{x) = Im \iiUp{x) where Up are the ordered 
products of the link variables Ui, around the pla- 
quette. In order to understand an eventual real- 
ization of the index theorem on the lattice one 
would like to establish a relation between the 
number of real eigenvalues (compare the discus- 
sion above) of the fermion matrix and the func- 
tional I'iU]. 

After inspecting spectra for several configura- 
tions one can conjecture the rule 

^ of real eigenvalues = 4 |:/| . (7) 

It has to be remarked, that this is not a result in a 
mathematically strict sense, since isolated gauge 
field configurations which violate (0) can be con- 
structed. However, it can be tested numerically 
whether (^ holds for a considerable portion of 
the gauge field configurations of the lattice path 
integral. In order to analyze this question we per- 
formed a simulation of lattice QED2 with two fia- 
vors of dynamical fermions at several values of /3 
(inverse coupling squared) and L(T = L). The 
value of K was always chosen close to the criti- 
cal K for the respective values of /3 and L. For 
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each Monte Carlo configuration the topological 
charge @ was evaluated. We define p(/3) to be 
the probability of finding (|^) correct for a given 
P (and L). 
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Figure 1. Probability p{f3) of finding (0) correct 
as a function of /?. We show our results for lat- 
tice sizes L ~ 4, 8, 12 and 16. The symbols are 
connected to guide the eye. 

It is obvious from fig. 1 that already for moder- 
ately high values of /? the configurations obeying 
entirely dominate the lattice path integral. 
This establishes, that with the proper interpre- 
tation of the real eigenvalues and their eigenvec- 
tors as the trace of configurations with non-trivial 
topology, the index theorem holds in a probabilis- 
tic sense. In particular in the continuum limit 
/? ^ oo (|7|) holds with probability 1. 

At this point we should point out the particular 
features of QED2 as opposed to QCD4. For QED2 
it is known Q], that in addition to the Atiyah 
Singer index theorem there holds the so-called 
vanishing theorem, which states that there are ei- 
ther only left-handed or only right-handed zero 
modes. Thus the result is exactly what one 
expects from the index and vanishing theorems 
for smooth continuum configurations (the factor 
4 comes from the extra real eigenvalues caused 
by the doublers). In QCD there is no vanishing 
theorem. Thus one expects the difference of the 
numbers of left- and right-handed real eigenvalues 
to be related to the topological charge. 



How can one utilize an understanding of the 
realization of the index theorem on the lattice? 
In general all semiclassical continuum arguments 
involving the index theorem can be extended to 
the fully quantized lattice model. In [Q it was 
e.g. demonstrated how the dependence of the 
pseudoscalar density on the topological charge 
can be understood using the index theorem. In 
addition one can use the results for the spectrum 
to analyze problems that are specific for the lat- 
tice approach. In particular in [Q we investigated 
the behaviour of the sign of the fermion deter- 
minant, establishing that it can become negative 
only for configurations with non-trivial topology. 

We believe, that extending the numerical part 
of our analysis to QCD4 will give interesting in- 
sight on the role of topologically non-trivial gauge 
field configurations in the fully quantized model. 
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